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1.  Introduction 

In  most  signal  processing  techniques,  the  elements  of  a  sensor  array  are  treated  as 
independent  entities.  However,  in  practice,  mutual  coupling  exists  between  the  array  sensors. 
Because  the  mutuals  change  the  sensor  impedances,  the  gain  and  radiation  pattern  of  the  array 
can  be  greatly  distorted.  Gupta  and  Ksienski  [1]  investigated  the  effects  of  mutual  coupling 
on  the  performance  of  an  adaptive  array.  In  their  treatment,  the  matrix  Z0  characterizing  the 

mutual  coupling  between  the  sensors  was  determined  using  a  mathematical  model  which 
models  the  antenna  array  consisting  of  J  sensors  as  an  (J  + 1)  terminal  network.  This  matrix 
was  then  used  to  determine  the  sensor  outputs  that  would  have  existed  had  there  been  no 
mutual  coupling.  A  compensation  scheme  was  then  developed  to  study  the  performance  of 
the  conventional  beamforming  technique.  Yeh  and  Leou  [2]  used  the  same  mathematical 
model  and  applied  it  to  the  MUSIC  algorithm.  Recall,  in  the  MUSIC  algorithm,  that  one 
plots  the  inverse  of  the  correlation  between  the  noise  subspace  Enand  the  directional  vector 

a(0) ,  which  we  denote  by|||Ena(0)||2  j  .  If  mutual  coupling  is  present,  Yeh  and  Leou  show 

.  Failure  to  do  so  results  in  severe 

degradation  of  the  estimates.  Shau  [3]  considered  the  case  of  deterministic  signals  in  a  noise 
free  environment  and  eliminated  the  effects  of  mutual  coupling  for  the  method  known  as  the 
Modified  Forward  Backward  Linear  Prediction  (MFBLP).  Himed  [4]  considered  the  case  of 
mutual  coupling  for  narrowband  signals  and  effectively  compensated  for  the  effects  of  mutual 
coupling  for  the  techniques  known  as  the  Matrix  Pencil  Approach,  with  two  different 
operators,  namely  the  Moving  Window  and  ESPRIT.  Two  recent  analyses  [6,7]  have  also 
considered  the  effects  of  mutual  coupling.  In  [6],  the  performance  of  high-resolution 
direction  finding  systems  such  as  MUSIC  is  studied  under  a  wideband  scenario.  They  show 
that  performance  is  adversely  affected  by  mutual  coupling  not  only  at  the  low  end  of  the 
frequency  band  but  also  at  the  high  end.  They  devised  ways  to  correct  the  actual  voltage 
matrix  using  the  terminal  impedance  matrix  that  is  derived  from  the  MOM  impedance 
matrix.  Their  results  approach  the  ideal  case.  However,  this  correction  does  not  involve  the 
alteration  of  the  signal-processing  algorithm.  In  [7],  an  adaptive  algorithm  that  accounts  for 
the  mutual  coupling  between  the  elements  of  an  array  was  studied.  It  is  shown  that  failure  to 
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account  for  mutual  coupling  causes  the  failure  of  all  adaptive  techniques.  Their  technique  is 
based  again  on  the  MOM  admittance  matrix  to  quantize  the  mutual  coupling  and  to  find  the 
optimum  weights  to  produce  the  desired  nulling. 

In  this  effort  we  develop  algorithms  to  effectively  compensate  for  the  effects  of 
mutual  coupling  when  using  space-time  adaptive  processing  techniques.  However,  for 
simplicity,  we  will  only  use  the  model  described  in  [4].  The  next  section  describes  the 
proposed  model. 

2.  Model 

We  restrict  ourselves  to  a  simple  model  and  consider  a  linear  array  of  J  dipoles 
uniformly  spaced  at  a  distance  d.  Each  dipole  is  of  length  t  and  has  a  radius  r  satisfying  the 
condition  r«t.  A  load  is  attached  to  the  center  gap  of  each  dipole.  Assume  there  are  q 
narrowband  signals  impinging  on  the  array  as  planar  wavefronts.  The  voltages  induced  by 
the  assumed  signals  on  the  loads  are  the  outputs  of  the  dipoles.  Induced  currents  will  appear 
on  the  dipoles.  These  currents  reradiate  and  generate  scattered  fields.  The  scattered  fields 
then  induce  currents  on  the  neighboring  dipoles.  The  process  of  induction  and  reradiation 
causes  the  mutual  coupling  among  the  dipoles. 


Using  one  sinusoidal  expansion  and  weighting  function  per  dipole,  the  method  of 
moments  [5]  was  used  to  obtain  the  matrix  of  mutuals  (see  Figure  1). 
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Figure  1.  Sub-sectioning  of  the  J  dipoles. 
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Denote  the  current  distribution  by  k(z)  ,  (assuming  longitudinal  distribution  and 
neglecting  all  other  distributions)  and  the  i-th  expansion  function  by  fj(z).  Then 

k(Z)=E!(*(z)  (1) 

i=l 

where  l(i) ;  i=l,2,...,J ,  denotes  the  unknown  current  amplitude  to  be  determined  on  each 
dipole.  At  a  point  (y,z)  in  the  Y  -  Z  plane,  the  scattered  field  is  given  by 

Ew(y,z)=i;i(iW')(y.z)  (2) 

i=l 

where  E(s)(y,z)  is  the  scattered  field  from  the  i-th  dipole.  The  total  field  will  then  be 


E(y,z)=E(inc)(y,z)+E(s)(y,z)  (3) 

where  E(inc)  is  the  incident  field.  Let  Ez  be  the  z-component  of  the  total  field.  A 
generalized  voltage  V(k)  induced  on  the  subsection  spanned  by  the  function  fk(z)  can  be 
defined  with  respect  to  a  weighting  function  wk(z)  as 

v(k)=F|EI(y,z),wk(z)]  (4) 

where  F(.,.)  is  bilinear  with  respect  to  E2(y,z)  and  wk(y,z).  Similarly,  we  define  the 
voltage  produced  by  the  incident  field  on  the  k-th  dipole  by 

V(inc)(k)  =  F[E(zi,,c)y,z),wk(z)]  (5) 

and  the  voltage  produced  by  the  scattered  field  on  the  k-th  dipole  by 

V(s)(k)  =  F[E[s)y,z),wk(z)].  (6) 

Thus,  the  total  voltage  introduced  in  the  k-th  dipole  is 

V(k)=V(inc)(k)+V(s)(k),  (7) 

which,  for  metallic  scatterers,  becomes 

V(k)  =  VM(k)+  V(s)(k)=  0  (8) 


or 


V(inc)(k)=  -  V(s)(k). 


(9) 


However, 
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(10) 


V(s)(k)  =  F  ^l(i)ESs)(y,z)wk(z)  =  £l(i)F[E[s)(y,z),  wk(z)]. 

_i=l  J  w 

The  total  impedance  between  the  k-th  and  i-th  dipoles  is  defined  to  be 

z*  =-F[E1(s)(y,z)wk(z)].  (11) 

Thus, 

V<s)(k)  =  £-  z“l(i)  ;  k  =  l,2,-,J.  (12) 

i=l 

In  matrix  notation 

V(s)=-ZI  (13) 

where 

V{,)  =  [v(s)(i)  ,  V(s)(2)  ,  ...  ,  V(s)(j)f  (14) 

and 

I  =  [l(l)  ,  1(2)  ,  ».  ,  I(J)]T.  (15) 

The  total  impedance  matrix  Z  can  be  decomposed  into  two  parts  as 

Z=Z0+ZL  (16) 

where,  Z0  is  the  generalized  impedance  matrix  and  ZL  is  the  load  matrix. 

Assuming  that  all  loads  are  loaded  with  the  same  load  z, ,  the  matrix  Zl  is  a  diagonal  matrix 
with  elements  z,  and  is  given  by 

ZL=diag(z,).  (17) 

The  ki-th  element  of  Z,  therefore,  is 

Z*”  =  Zy  +  Z|8y  ,  (18) 

where  Zy  is  the  mutual  impedance  between  the  k-th  and  i-th  dipoles.  The  total  voltages 
induced  on  a  load  z,  are  given  by 

V(l)  =  ZlI  and  I=Z‘'V(L).  (19) 

However, 

v(inc>  =ZI=Z0Z[' V(L)  +V(L)  =(i  +  Z0Zl  )v(L),  (20) 

which  implies  that 
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V(L)  =(l+Z0ZL1)'IV(inc).  (21) 

Let  H  be  the  matrix 

H  =  (i+Z0Zl).  (22) 

H  can  be  written  as 

l  +  (zn/zi)  (z12/z,)  •••  (zu/z,) 

(z2i/z,)  l  +  fzjj/z,)  -  (z2J/z,) 

(ZJl/Zl)  (Z  J2  /Z1 )  l  +  (ZJj/Zl 

Thus,  when  incident  signals  are  impinging  on  the  array  and  in  the  presence  of  additive  noise, 
the  output  of  the  linear  array,  at  time  n,  will  be 

V(L)(n)  =  H"1  V(iDc)(n)  +  w(n) ,  n=l,2,  ...,N  (24) 

where  N  is  the  number  of  pulses  and  w(n)  is  the  additive  white  noise  vector. 

Note  that  usually,  the  vector  V^(n)  consists  of  signal  vectors,  clutter  and  jamming  signals. 
The  problem,  therefore,  becomes  one  of  accounting  for  mutual  coupling  in  the  data 
generation  and  ways  to  compensate  for  these  effects  in  the  signal  processing  stage. 

Of  particular  interest  to  us  is  a  signal  generation  code  developed  by  Scientific  Studies 
Corporation  (SSC)  for  airborne  radar  applications  [8].  This  data  generation  technique  is 
being  used  in  the  Air  Force  Research  Laboratory  Multi-Channel  Signal  Processing  System 
(AFRL-MCSPS).  In  the  remainder  of  this  report,  we  will  refer  to  this  section  of  the  code  as 
the  physical  model. 

The  physical  model  is  a  very  powerful  tool  for  signal  generation,  especially  for 
algorithm  evaluation  and  comparison.  It  provides  the  user  with  both  a  datacube  consisting  of 
J  channels,  N  pulses  and  K  range  cells,  and  the  true  covariance  matrices  of  all  components 
involved  in  the  simulation  (target,  clutter  and  jamming).  This  allows  for  optimal 
performance  evaluation  of  STAP  algorithms.  However,  for  our  purposes,  we  modified  the 
code  so  as  to  generate  the  space-time  true  covariance  matrices  as  well  as  estimates  of  these 
matrices  generated  using  the  sample  covariance  matrix. 
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3.  Modification  of  the  Physical  Model 

The  first  task  in  this  effort  is  to  incorporate  the  effects  of  mutual  coupling  into  the 
MCSPS’s  physical  model  (data  generation).  Therefore,  care  must  be  taken  to  incorporate 
these  effects  into  all  of  the  signals  involved  in  the  simulation,  i.e.;  target,  clutter  and 
jamming. 

Note  that  the  physical  model  [8]  enables  the  user  to  generate  a  sequence  of  data  as 
well  as  the  “true  covariance”  matrices  of  the  involved  signals  [8].  First,  we  made  changes  to 
the  code  so  that  the  output  sequence  will  be  available  to  the  user  with  each  component 
identified  separately.  This  is  desired  for  diagnostics  purposes. 

A  MATLAB  code  was  written  to  compute  the  matrix  of  mutuals.  This  was  done  for 
the  case  of  a  linear  uniformly  spaced  array  composed  of  have-wavelength  dipoles  as 
described  in  section  2.  The  load  impedance  zL  was  assumed  to  be  z j, .  This  then  was 
incorporated  into  the  physical  model  and  the  user  can  now  generate  space-time  data  with  the 
effects  of  mutual  coupling  included  in  the  data  generation.  At  this  stage,  the  matrix  of 
mutuals  is  applied  to  all  data  vectors,  however,  it  is  also  possible  for  the  user  to  choose  which 
data  sets  (i.e.;  targets,  jammers  or  clutter)  the  matrix  of  mutuals  will  be  applied  to. 

To  illustrate  this,  consider  the  following  signal  vector  s,  usually  given  by 

s=a(st®s5),  (25) 

where  a  is  a  constant  (complex)  amplitude,  s,  is  a  temporal  steering  vector,  ss  is  a  spatial 

steering  vector  and  ®  denotes  the  Kronecker  tensor  product.  In  the  presence  of  mutual 
coupling,  only  the  spatial  steering  vector  is  affected  and  is  then  modified  as  follows 

Sm-H-'s..  (26) 

The  new  signal  vector  sz  then  becomes 

sz=ai®s5,z)=a|st®(H-,sJ.  (27) 

A  similar  operation  is  applied  to  all  signal  vectors  including  target,  clutter  and  jamming. 

This  does  not  affect  additive  white  noise. 

Consider  now  the  received  space-time  vector,  which  we  denote  by  X-  X  can  be 
written  as 
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XT  =  [x(l)  x(2)  -  x(n)],  (28) 

where  the  vector  x(n)  ;  n=l,2,...,N ,  is  a  composite  of  target  signal,  clutter,  jamming  and 
additive  white  noise.  It  is  generally  given  as 

x(n)=t(n)+c(n)+i(n)+w(n)  ;  n  =  1, 2, ...,  N  (29) 

where  t(n)  is  the  target  signal  vector,  c(n)  is  the  clutter  signal  vector,  i(n)is  the  jamming 
signal  vector  and  the  w(n)is  the  additive  white  Gaussian  noise  assumed  to  be  independent  and 
identically  distributed. 

In  the  absence  of  mutual  coupling,  the  covariance  matrix  of  the  received  signal  vector 
X  is  given  by 

Rx  =E{xXH}=Rt+Rc+Ri+RD,  (30) 

where  Rt ,  Rc ,  Rj  and  Rn  are  the  covariance  matrices  of  the  target,  clutter,  jamming  and 
noise,  respectively. 

To  further  simplify  our  analysis,  we  assume  that  the  datacube  consists  of  J  channels, 

N  pulses  and  K  range  cells.  For  simplicity,  let  y  denote  a  general  vector  used  to  denote  any 

specific  process,  t,  c,  or  i.  The  received  space-time  vector  y  has  therefore  dimension  JNxK. 
Each  column  of  y  is  assumed  to  have  been  formed  by  concatenating  all  channels  for  pulse  1 
through  pulse  N.  Let  y(n)  be  the  Jxl  signal  vector  obtained  at  pulse  n. 

yTW  =  [yi(n)  y2(n)  •••  y,(n)],  (31) 

where  yk(n)  is  the  signal  received  at  the  k-th  channel  and  n-th  pulse.  Then,  the  space-time 
vector  Y  can  be  expressed  as  (in  concatenated  form) 

YT=k(l)  y(2)  ...  y(N)].  (32) 

In  the  presence  of  mutual  coupling,  the  vector  y^(n)  replaces  the  vector  y(n) ,  where 
yz(n)  is  given  by 

yz(n)=H_1y(n)  ;  n  =  1, 2, ....,  N.  (33) 

The  received  space-time  vector  is  then  changed  to 

Yl={yz(l)  yfi)  ...  yz(N)}=  (H-'y(l)  H-'y(2)  ...  H-1y(N)}.  (34) 
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Let  T  be  the  following  matrix 


H'1  0  •••  0 

0  H"1  •••  0 

0  0  •••  H_1 

Then  Yz  can  be  expressed  as 

YZ=T Y. 


r=diag(H-')= 


(35) 


(36) 


Hence,  in  the  presence  of  mutual  coupling,  the  received  signal  vector  is  expressed  as 
xx(n)=H_1[t(n)+c(n)+i(n)]+w(n)  ;  n=  1,2,  ...,N, 
and  the  concatenated  space-time  vector  can  be  written  as 

xz=rx+w,  (37) 

where  W  is  the  space-time  additive  white  noise  vector.  Thus,  the  received  signal  (target, 
clutter,  jamming  +  noise)  covariance  matrix  is  given  by 

R...  -  }-  r(R,  +RC  +Ri)r"  +  Rw .  (38) 


Therefore,  in  the  presence  of  mutual  coupling,  Equations  (37)  and  (38)  are  used  throughout 
the  simulation  code. 

For  diagnostics  purposes,  we  applied  the  matrix  of  mutuals  to  each  component  (target, 
jamming  and  clutter)  individually  and  made  all  signal  components  available  to  the  user.  We 
have  also  included  the  effects  of  the  mutual  coupling  on  the  true  covariance  matrices  of  the 
targets,  jamming  and  clutter,  as  descried  by  Equation  (38). 

We  now  illustrate  the  impact  of  mutual  coupling  upon  both  the  true  covariance  matrix 
as  well  as  the  sample  matrix  estimator  defined  as 

•  (39) 

K-  k=l 


where  Xk  is  the  space-time  composite  signal  (target,  clutter,  jamming  and  noise)  vector.  Let 
R  be  the  true  covariance  matrix. 

We  consider  a  simulation  for  a  linear  array  consisting  of  J  =  14  channels  uniformly 
spaced  at  d  =  (^/2)  and  a  coherent  processing  interval  (CPI)  of  N  =  18  pulses.  A  pulse 
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repetition  frequency  (PRF)  of  300  Hz  is  assumed  and  a  450  MHz  transmit  frequency  is  used. 
In  the  scenario,  we  had  a  target,  two  jammers,  clutter  and  additive  white  noise. 

Figures  2  and  3  show  the  output  of  the  real  and  imaginary  part,  respectively,  of  the 
true  channel  1  temporal  covariance  sequence  without  and  with  mutual  coupling  as  a  function 
of  the  temporal  lag.  This  sequence  is  obtained  from  the  (1,1)  element  of  each  JxJ  block 
matrix  contained  in  the  first  block  row  of  the  JNxJN  covariance  matrix. 


Figure  2.  Real  part  of  true  temporal  covariance  sequence  for  channel  1 . 
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Figure  3.  imaginary  part  of  true  temporal  covariance  sequence  for  channel  1. 


From  Figures  2  and  3,  note  that  mutual  coupling  does  change  the  true  temporal 
covariance  sequence.  The  change  is  more  accentuated  in  the  imaginary  part.  Even  though, 
mutual  coupling  affects  the  spatial  component  of  the  covariance  matrix,  the  pre-  and  post¬ 
multiplications  by  the  matrix  of  mutuals  of  the  space-time  covariance  matrix  do  affect  the 
temporal  aspect  of  the  matrix. 

Figures  4  and  5  show  the  output  of  the  real  and  imaginary  part  of  the  sample 
covariance  matrix,  estimated  from  448  secondary  data  cells. 
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Figure  5.  imaginary  part  of  temporal  sample  covariance  sequence. 


As  noted  in  Figures  4  and  5,  mutual  coupling  also  affects  the  sample  covariance 
matrix  estimator.  This  is  expected  since  the  matrix  of  mutuals  affects  the  sample  covariance 
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matrix  in  the  same  manner  as  the  true  covariance  matrix.  The  differences  between  the  true 
and  estimated  values  are  due  to  the  small  sample  support  size  (K)  used  in  the  estimation 
procedure  of  Equation  (39). 

Figures  6  and  7  show  the  spatial  variations  of  the  real  and  imaginary  parts, 
respectively,  of  the  true  space-time  covariance  matrix  with  and  without  mutual  coupling.  In 
these  figures,  we  plot  the  cross-channel  correlation  sequence  as  a  function  of  the  spatial  lag 
(i.e.,  the  integer  difference  between  the  channel  numbers)  with  channel  1  used  as  the 
reference  channel).  We  emphasize,  however,  that  spatial  stationarity  does  not  hold  in  the 
presence  of  mutual  coupling  between  elements.  Therefore,  we  would  expect  variations  in  the 
cross-channel  correlation  sequence  when  using  another  channel  as  the  reference  channel. 
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Figure  7.  Imaginary  part  of  true  spatial  covariance  sequence. 


Figures  6  and  7  are  important  because  they  clearly  show  the  effects  of  mutual  coupling.  As 
we  have  stated  earlier,  mutual  coupling  affects  the  spatial  variation  aspect  of  both  the 
received  signal  vector  as  well  as  its  hue  covariance  matrix.  This  shows  very  clearly  in  both 
figures  as  the  real  and  imaginary  components  of  the  hue  covariance  matrix  deviate  from  their 
values  when  mutual  coupling  is  included. 

Figures  8  and  9  show  the  real  and  imaginary  parts  of  the  estimated  sample  covariance 
matrix.  The  procedure  used  here  is  described  in  the  MATLAB  code  given  above.  For 
simplicity,  only  the  first  J  terms  of  the  first  column  are  used  in  the  simulation. 
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0  2  4  6  8  10  12  14 

Spatial  Lag,  I 

Figure  8.  Real  part  of  spatial  sample  covariance  sequence. 


Figure  9.  Imaginary  part  of  spatial  sample  covariance  sequence. 


Again,  the  same  effects  are  observed  in  this  case. 
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Power  Spectrum,  S(fd,fs)  (dB) 


The  next  set  of  figures  shows  both  the  true  power  spectrum  and  that  obtained  through 
the  estimated  sample  covariance  matrix,  without  mutual  coupling  first  and  then  with  mutual 
coupling  included. 


Without  mutual  coupling _ With  mutual  couplin 


Figure  10.  Normalized  true  power  spectrum. 


Figure  11.  Normalized  estimated  power  spectrum. 
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Figure  12.  Normalized  true  power  spectrum  (modified  periodogram). 


Note  from  Figures  10-12  that  mutual  coupling  does  change  the  respective  power 
spectra.  In  particular,  it  can  be  seen  from  figure  10  that  the  angle-Doppler  map  has 
noticeably  changed  when  mutual  coupling  is  present.  Mutual  coupling  introduces  a  strong 
ridge  across  the  spatial  frequency  domain  at  zero  Doppler.  Note  also  from  figure  1 1  that  a 
strong  ridge  appears  at  zero  Doppler.  The  same  effect  is  also  seen  in  figure  12. 

4.  Application  to  Space  Time  Adaptive  Processing 

In  section  3,  we  included  the  effects  of  mutual  coupling  into  the  physical  model  of  the 
MCSPS,  from  which  data  is  generated  for  future  processing.  We  have  shown  how  mutual 
coupling  changes  the  received  signals  as  well  as  their  respective  covariance  matrices.  In  this 
section,  we  apply  some  well-known  STAP  algorithms  such  as  the  Adaptive  Matched  Filter 
(AMF),  also  referred  to  as  the  CFAR-AMF,  to  determine  the  effects  of  mutual  coupling  on 
the  performance  of  these  algorithms.  To  achieve  this  goal,  we  generate  two  sets  of  data. 

First,  only  the  true  space-time  covariance  matrices  are  generated.  These  in  turn  will  be  used 
in  the  signal  processing  aspect  of  the  algorithm.  Then,  a  datacube  consisting  of  J  channels,  N 
pulses  and  K  range  cells  is  generated.  The  needed  covariance  matrices  are  estimated  using 
this  data.  We  first  apply  the  CFAR-AMF  without  mutual  coupling  and  determine  the 
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performance  of  the  algorithm  in  these  conditions.  The  metrics  used  to  determine  the 
performance  of  the  algorithm  are  the  output  SINR  and  the  adapted  antenna  beam  pattern, 
both  plotted  as  a  function  of  angle  and  Doppler. 

The  physical  model  was  modified  so  that  true  space-time  covariance  matrices  and  a 
datacube  are  generated.  The  covariance  matrices  have  dimension  (JNxJN)  and  the  datacube 
has  dimension  (JNxK).  We  have  successfully  been  able  to  generate  “true”  and  “estimated” 
covariance  matrices  suitable  for  STAP  processing.  This  is  very  practical  when  performing 
detection  analyses.  We  programmed  the  CFAR-AMF  approach  so  that  mutual  coupling 
effects  will  be  analyzed.  As  mentioned  earlier,  two  measures  have  been  used  to  study  the 
performance  of  the  algorithm,  which  include  the  estimated  Signal  to  Interference  plus  Noise 
Ratio  (SINR)  and  the  adapted  beam  pattern  plotted  in  the  angle-Doppler  domain.  Both 
measures  were  plotted  as  a  function  of  the  steering  vectors’  azimuth  and  Doppler.  The 
estimated  SINR  tells  us  how  well  we  have  recovered  the  target  after  cancellation  of  the 
interference,  which  in  turn  is  best  shown  by  the  adapted  beam  pattern.  These  are  defined  as 


SINRest  = 


vr-'xh 


VR-'Vh 


and 


Patadp=VR-'VH) 


(40) 


(41) 


respectively.  Note  that  the  adapted  pattern  is  just  the  denominator  of  Equation  (40).  We 
have  also  plotted  the  variations  of  the  numerator  of  Equation  (40),  sometimes  referred  to  as 
the  sample  matrix  inversion  (SMI)  test  statistic,  over  angle  and  Doppler.  It  is  given  by 


SMI  = 


vr-'xh 


(42) 


The  scenario,  which  we  ran,  was  similar  to  the  mountaintop  data  set  [9].  It  consists  of  J=14 
channels,  N=16  pulses  and  K  =  3(JN)=672  range  cells.  These  will  be  useful  for  covariance 
computation.  One  target  with  a  signal  to  noise  ratio  (SNR)  of  -5  dB,  361  clutter  patches 
constituting  a  clutter  to  noise  ratio  (CNR)  of  38  dB,  2  jammers  with  a  jammer  to  noise  ratio 
(JNR)  of  47  dB  have  been  included  in  the  simulation.  The  target  has  a  velocity  of  33.3  m/s 
and  is  located  at  an  azimuthal  angle  of  0°  (mainbeam).  This  translates  into  a  normalized 
Doppler  frequency  of  0.333  and  a  normalized  spatial  frequency  of  zero.  Two  “true” 
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covariance  matrices,  one  without  mutual  coupling  and  the  other  with  mutual  coupling  were 
generated.  We  also  generated  two  data  vectors  at  the  cell  under  test,  one  without  mutual 
coupling  and  the  other  with  mutual  coupling.  Let  R,  V  and  X  denote  the  interference 
covariance  matrix,  the  steering  vector  and  the  data  vector,  respectively,  in  the  absence  of 
mutual  coupling.  In  the  presence  of  mutual  coupling,  these  quantities  are  denoted  by  Rz , 
Vz  and  Xz ,  respectively,  where  Vz  is  defined  below  in  Equation  (43).  We  have  generated 
three  sets  of  plots.  In  the  first  case,  no  mutual  coupling  was  included.  In  the  second  case, 
mutual  coupling  is  included,  however  the  steering  vector  in  Equations  (40)-(42)  did  not 
change.  In  the  third  case,  mutual  coupling  was  included  and  the  steering  vector  has  been 
replaced  by  Vz .  This  is  shown  below. 

•  Case  LSINR^ 

•  Case  2:811^, 

•  Case  3:  SINR^, 


VR-‘Xh 

2 

2 

vr‘vh 

-,  Patadp  =  VR-'V11  and  SMI= 

VR’X 

vr;'x“ 


vr;‘v 


r-J-,  Patadp  =  VR;‘Vh  and  SMI  =  VR;'X 


Yzr;‘Xz 


ij-,  Pat.*,  =  Y.RI'V"  and  SMI  =  |V,R;‘X 


V.RZ’V 


I  -1 ' 


2  Z  -I_z 
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Figure  14.  Adapted  beam  pattern,  (a)  without  mutual  coupling,  (b)  with  mutual  coupling,  and 

(c)  with  compensating  for  mutual  coupling. 
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Figure  15.  Test  statistic,  (a)  without  mutual  coupling,  (b)  with  mutual  coupling,  and  (c)  with 

compensating  for  mutual  coupling. 


In  the  ideal  case,  (Figures  13-15,  case  (a)),  where  no  mutual  coupling  was  assumed, 
the  CFAR-AMF  processor  clearly  distinguishes  the  target  and  cancels  the  interference.  A 
sharp  peak  is  generated  at  the  location  of  the  target,  and  two  nulls  have  been  generated  at  the 
locations  of  the  jammer  and  a  null  is  generated  along  the  clutter  ridge.  When  mutual 
coupling  was  involved,  and  using  ideal  steering  vectors,  (Figures  13-15,  case  (b)),  the  target 
is  completely  buried.  It  is  extremely  difficult  to  identify  the  target,  even  though  appropriate 
nulls  have  been  generated  at  the  jammer  locations  and  along  the  clutter  ridge.  The  only 
difference  is  that  the  nulls  are  broader  and  less  deep.  This  clearly  shows  that  mutual 
coupling  degrades  the  performance  of  the  CFAR-AMF.  Therefore,  appropriate  measures 
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have  to  be  developed  so  as  to  compensate  for  these  effects.  At  this  stage,  our  compensating 
scheme  consists  of  applying  steering  vectors  corrupted  with  mutual  coupling.  As  discussed 
earlier,  the  new  steering  vector  has  the  form 

Yz=rV,  (43) 

where  V  is  the  ideal  steering  vector  and  T  is  the  following  matrix 


r= 


H-1 

0 


0 

H'1 


0  0 


0 

0 


(44) 


Note  from  Figures  (13)-(15),  case  (c),  that  this  simple  compensating  scheme  worked  well  in 
the  sense  that  the  target  has  been  recovered  better  and  that  the  nulls  are  also  clearly 
distinguished. 

The  intent  of  this  report  was  to  include  the  effects  of  mutual  coupling  into  our 
simulation  capabilities  for  future  analyses  regarding  their  effects  on  STAP  techniques. 

Future  research  should  be  devoted  to  the  study  of  other  sub-optimum  techniques.  These 
would  include  both  factored  approaches  (FTS  and  FST)  [10],  the  extended  factored  approach 
(EFA)  [11],  adaptive  displaced  phase  center  aperture  (ADPCA)  [12],  joint  domain  localized 
(JDL)  [13],  the  principal  components  technique  (PC)  [14],  the  cross-spectral  metric  (CSM) 
[15]  and  the  newly  developed  technique  termed  as  the  parametric  adaptive  matched  filter 
(PAMF)  [16]. 


5.  Effects  of  Mutual  Coupling  on  the  True  SINR 


We  have  also  investigated  the  performance  of  the  sample  covariance  matrix  estimator 
and  compared  the  performance  of  the  output  SINR  as  a  function  of  number  of  the  sample 
support  size.  The  output  SINR  is  given  by: 


SINR= 


vhr_ix 

2 

" 

VhR_1RR~'V  ’ 


(45) 
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A 

where  R  is  the  true  covariance  matrix  and  R  is  the  sample  covariance  matrix  estimated  using 
K=nJN  secondary  data.  Note  that  n  is  an  integer,  J  is  the  number  of  channels  and  N  is  the 

number  of  pulses. 

In  the  absence  of  mutual  coupling,  the  following  table  shows  the  difference  between 
the  true  SINR  and  the  estimated  SINR  as  a  function  of  n. 


Table  1.  SINR  as  a  function  of  n,  no  mutual  coupling 


n 

20 

60 

SINR,  (dB) 

28.3572 

28.3572 

SINRe  (dB) 

26.5697 

26.5889 

In  the  presence  of  mutual  coupling,  the  following  table  shows  the  difference  between 
the  true  SINR  and  the  estimated  SINR  as  a  function  of  n. 

Table  2.  SINR  as  a  function  of  n,  with  mutual  coupling 


1  n 

20 

60 

1  SINR,  (dB) 

24.2185 

24.2185 

I  SINRe  (dB) 

22.4393 

22.4457 

The  above  tables  demonstrate  two  things: 

1 .  There  is  a  4  dB  difference  between  the  coupled  and  uncoupled  SINR,  for  both  the 
true  and  estimated  SINR. 

2.  The  estimated  SINR  approaches  the  true  SINR  as  n  increases. 

6.  Probability  of  Detection 

In  this  section,  we  investigate  the  effects  of  mutual  coupling  on  the  probability  of 
detection  (Pd)  of  the  adaptive  CFAR-AMF  algorithm  for  a  pre-determined  probability  of 
false  alarm  (Pfa).  A  single  target  with  an  SNR  of  -5dB  was  assumed.  We  first  checked  some 
prior  results  from  an  AFRL  in-house  analysis  conducted  by  Dr.  James  Michels.  After 
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conducting  validation  analyses,  the  MATLAB  codes  developed  here  were  optimized  to 
reduce  computation  run  times  and  were  then  used  to  generate  probabilities  of  detection  (Pd) 
curves  for  several  cases,  using  Monte-Carlo  procedures.  We  noticed  significant  speedup 
improvements  of  the  MATLAB  code  that  we  implemented.  Next,  we  included  the  effects  of 
mutual  coupling  and  then  computed  the  corresponding  probability  of  detection. 

Figures  ( 1 6)-(2 1 )  show  the  performance  of  the  CFAR-AMF  without  and  with  mutual 
coupling  under  different  scenarios.  The  probability  of  false  alarm  is  kept  unchanged  for  all 
cases  under  study  and  is  equal  to  0.01  (Pfa  =  0.01).  For  simplicity  and  speed  up,  in  the  first 
case,  J  =  2  channels  and  N  =  2  pulses  have  been  considered.  The  results  shown  below  have 
been  generated  using  10,000  Monte-Carlo  snapshots  and  averaged  over  1000  independent 
runs. 

Figure  16  shows  the  detection  performance  of  the  CFAR-AMF  with  K  =  1000  range 
cells  used  to  estimate  the  disturbance  covariance  matrix.  The  normalized  Doppler  frequency 
is  assumed  to  be  zero  while  the  elevation  and  azimuth  angles,  <j)  and  0,  are  also  assumed  to  be 
zero. 


Figure  16.  Pd  vs.  SNR,  K  =  1000. 
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As  can  be  seen  from  Figure  16,  to  achieve  the  same  Pd  with  and  without  mutual  coupling,  an 
increase  of  about  5  dB  in  the  signal  to  noise  ratio  is  needed.  This  also  can  be  interpreted  as 
having  a  loss  of  about  4-5  dB  in  Pd  when  mutual  coupling  is  included  in  the  simulation. 

Figure  17  is  similar  to  Figure  16.  In  this  case,  however,  only  8  secondary  data  vectors 
have  been  used  for  covariance  matrix  estimation. 


CFAR-AMF 
Pd  vs.  SNR 


Figure  17.  Pd  vs.  SNR,  K  =  8. 


Figure  17  shows  very  similar  Pd  behavior  as  those  of  figure  16.  The  only  difference 
is  that  with  K  =  8,  the  covariance  matrix  is  estimated  with  only  2JN  samples.  This  is  usually 
the  minimum  number  of  samples  used.  As  expected,  the  performance  is  lower  than  with 
K=1000,  but  the  loss  in  Pd  for  both  cases  is  the  same;  i.e.;  when  mutual  coupling  is  included, 
there  is  a  4-5  dB  drop  in  SNR  for  achieving  the  same  Pd. 

Figure  18  shows  the  variation  of  Pd  as  a  function  of  azimuth  angle  0  of  the  target, 
while  keeping  the  SNR,  the  normalized  Doppler  frequency  and  the  azimuth  angle  constant. 
The  Probability  of  false  alarm  (Pfa)  is  still  kept  fixed  at  0.01  and  K  =  200  samples  are  used 
for  covariance  matrix  estimation.  The  target  was  assumed  to  have  an  SNR  of  10  dB. 
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CFARAMF 

Pd  vs  e 


Note  from  the  above  figure  that  Pd  drops  by  more  than  50%  when  mutual  coupling  is 
included  in  the  simulation. 

Figure  19  shows  the  variation  of  Pd  vs.  the  target’s  normalized  frequency.  K  =  200 
samples  are  used  for  covariance  matrix  estimation,  Pfa  was  kept  at  0.01  and  an  SNR  of  10  dB 
was  assumed.  In  this  case  0  =  0°. 
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Figure  19.  Pd  vs.  Normalized  Doppler  Frequency. 


Note  from  the  above  figure  that  in  the  presence  of  mutual  coupling,  the  probability  of 
detection  drops  by  more  than  50%.  This  is  a  substantial  drop  in  the  performance  that  should 
be  accounted  for  when  developing  schemes  to  counter  mutual  coupling. 

Figure  20  is  similar  to  Figure  19.  The  only  difference  is  that  0  =  30°. 
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Figure  20.  Pd  vs.  Normalized  Doppler  Frequency. 


The  probability  of  detection  drops  even  further  when  mutual  coupling  is  present.  With 
mutual  coupling  present,  Pd  is  constant  and  is  close  to  0.36  while  it  was  close  to  0.43  when  0 
was  equal  to  0°. 

Figure  21  is  again  similar  to  Figures  19  and  20.  In  this  case  0  =  60°. 
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In  this  case  the  probability  of  detection  is  close  to  0.29  when  mutual  coupling  is  present. 

Figure  22  show  the  performance  of  the  CFAR-AMF  without  and  with  mutual 
coupling  in  the  case  where  J  =  4  channels  and  N  =  8  pulses.  The  number  of  secondary  data 
used  is  K  =  64  (2JN  criterion).  The  probability  of  false  alarm  is  still  kept  unchanged  at  0.01 
(Pfa  =  0.01).  The  results  shown  below  have  been  generated  using  10,000  snapshots  and 
averaged  over  1000  runs.  The  normalized  Doppler  frequency  is  assumed  to  be  zero  while  the 
elevation  and  azimuth  angles,  0  and  <f>,  are  assumed  to  be  zero. 
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CFAR-AMF 

Pd  vs.  SNR 


Figure  22.  Pd  vs.  SNR,  K  =  64. 

As  can  be  seen  from  the  above  figure,  to  achieve  the  same  Pd  with  and  without  mutual 
coupling,  an  increase  of  about  5  dB  in  the  signal  to  noise  ratio  is  needed.  Therefore,  even 
when  J  and  N  are  increased  to  4  and  8,  respectively,  the  same  effect  is  seen. 

7.  Conclusions  and  Recommendations  for  Future  Work 

From  the  above  discussions  and  analyses,  it  can  be  concluded  that  mutual  coupling  is 
very  important  in  determining  the  performance  of  the  AMF.  It  can  be  seen  that  these  effects 
are  also  geometry  depended.  We  showed  that  Pd  drops  by  more  than  50%  depending  on  the 
case  studied.  In  some  other  cases,  we  showed  that  in  order  to  achieve  the  same  probability  of 
detection,  an  increase  of  about  4  dB  in  SNR  is  required.  It  is  our  recommendation  that  close 
attention  should  be  given  to  mutual  coupling  when  dealing  with  space-time  adaptive 
processors.  Care  should  be  taken  in  developing  compensation  methods  so  that  performance 
loss  is  minimized. 
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